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ABSTRACT 

The elastic transient response of an impulsively loaded circular 
ring is examined. The numerical stability criterion associated with 
a finite difference formulation of the problem is briefly discussed, 
and a method to reduce the stability restriction is presented. This 
method assumes that the density of the ring is greater in the tangen- 
tial direction than it is in the radial direction. The analytical 
solution to the ring equations with the artificial tangential inertia 
is determined for a speedfic impulsive load. This solution is 
compared with the solution to the ring equations with the correct 
tangential inertia. The comparison shows that on the basis of pre- 
dicting peak values of stress, the results for the ring with the 
artificial tangential inertia are very close to those of the unmod- 
ified ring. Thus, the method reduces the numerical stability require- 


ments and does not significantly affect the principal results. 
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CHAPTER I 


INTRODUCTION 


The finite difference method of analysis has long been : 
recognized as an effective and practical means for solving partial 
differential equations on a digital computer. The accuracy asso- 
ciated with this method can be adjusted by changing the incremental 
sizes of the independent variables involved. In the limit, when 
the sizes of the increments approach zero, the solution to the 
finite difference equations becomes identical to the solution to 
the corresponding differential equations. Consequently, the engineer 
may improve the accuracy of his finite difference solution by choos- 
ing sufficiently small values for the incremental sizes of the inde- 
pendent variables, with greater accuracy accompanied by increased 
execution time on the digital computer. However, there are many 
problems where the finite difference method imposes a restriction 
on the sizes of the increments. This restriction is known as the 
stability criterion and arises from a purely numerical phenomenon. 
The stability criterion defines a critical limit for the ratio of 
the increment sizes above which a numerically stable, valid solu- 
tion will result. A value of the ratio below this limit will lead 
to a divergent, meaningless solution. 

The stability criterion for the finite difference dynamic 


analysis of the thin, homogeneous, elastic, circular ring has been 


ly 


defined by sae Tt can be simplified to 


AA Da 2 
Tir > 2 Ga < | 





where A is thle radius of the rimg; h is the thickness, ofthe ring; 

D = Eh, the stretch rigidity; a is the ring density; and, Ad 
and f$\7T are the angular and time increments, respectively. A 
meaningful numerical solution for the transient response of the ring 
is possible only when the above restriction on A®D and AT 
is observed. A physical interpretation of the stability criterion 
says that the finite difference solution must propagate around the 
ring at least as rapidly as a stress wave. 

Once the size of Ad is chosen, the maximum size for the 


time increment becomes 


- (AP) A Ch (1) 


(AT) D 


MAX 


An accurate solution for the ring response requires a small value 
for the angular increment, AD . Therefore, according to Equation 
(1), eal eee is automatically restricted to a small value; and 
hence, large computation times will be required on the computer. 

To reduce the restrictions on the numerical stability of the 
ring response a method is proposed whereby the maximum time increment 


may be increased without producing an unstable solution. According 


WA 


Foubcuat tonal) «qm uss Jmay is determined once Ad is chosen. 
However, if the density -P is considered to be a variable in the 
stabiliiey critérion then™@ A, AAR will be proportional to the 
square root of the density at a particular value for Ad - The 
density in Equation (1) is a measure of the ring's inertia in the 
direction of the propagation of stress disturbances around the ring, 
that is, the tangential direction. The method proposed here modifies 
the ring by introducing an artificial tangential density that is 
larger than the actual ring density. Thus, the ring inertia will be 
characterized by the density ia in the radial direction and by a 
different density, ve » in the tangential direction. At» a given 
value of A a0) , the maximum time increment for a stable response 
of the modified ring will be proportional to the square root of 5P 
and, hence, will be larger than the critical time increment for the 
original ring by the factor vy. 

For example, for § = 4 the time increment for stable response 
of the modified ring is twice as large as the time increment for the 
original ring for the same Ag - Thus, the computation time can 
be cut in half. The decreased computational time afforded by the 
proposed method is a significant improvement on the conventional 
finite difference method of analysis. However, the introduction 
of the artificial tangential inertia will also lead to an incorrect 
solution for the ring response. Thus, the value of the method will 


depend upon the accuracy of the solution of the modified ring. 


lig 


This paper investigates the effect of an artificial tangential 
inertia on the analytical solution for a particular ring problem. 
Solutious are obtained for the elastic response of an impulsively 
loaded circular ring with and without artificial tangential inertia. 
The maximum stress in the ring is computed for the original ring 
( ¥ = 1) and for rings with ¥Y * 2, 3, and 4. A comparison of the 
maximum stress for each ring demonstrates the value of the method. . 

This study was completed during the 1967-1968 academic year at 
the Naval Postgraduate School, Monterey, California. Acknowledge- 
ment is gratefully made to Professor Robert E. Ball of the Aeronautics 


Department for his guidance and consultation. 
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CHAPTER II 


A DISCUSSION OF THE RING EQUATIONS 


This thesis considers a thin, circular ring composed of a 
homogeneous material. The governing equations are derived in 
Reference 4 and are repeated here for completeness. Figures 1 and 
2 illustrate the general geometrical characteristics. The radius A 
is measured to the middle surface of the ring, which lies at a dis- 
tance £ from both edges. A Sete ronbta element On tne tines d@ 
is located at an angle b from the bottom of the ring, where d is 
measured positive counter-clockwise. The displacement of a point 
on the middle surface is described by a radial and a tangential 
component, W and V, respectively. The radial distance from the 
middle surface to some arbitrary point on the ring is given by Z in 
Figure 2. 

A two-dimensional problem will be assumed with regard to the 
internal stresses of the ring. Figure 3 illustrates the forces 
acting on the element dad due to the normal and shear stress dis- 
tributions on the cross-sections. The forces N, M, and Q are, 
respectively, the normal force, the bending moment, and the radial 
shear force. These forces are stress-resultants referred to the 
middle surface of the element. The definitions of N and M are given 


by the following expressions 


h 
= 2 
ws oe Zz 
N = Oy dz = E é,dzZ ed 
Bn 4 


leo 


aw 
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GENERAL RING GEOMETRY 
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TWO-DIMENSIONAL STRESS-RESULTANTS 


(2b) 


s Zaz: [ees BER 
an 
2 


where Oy and Ce are the stress and strain on the cross section 


at a distance Z fromeene middie surface, 


In order to obtain an expression for the strain on in terms 


of Z, two additional assumptions are made: 


A line normal t6O™themmadd lewesummaice BDefore deformation 


Ls 
will be normal to this surface after deformation. 

2. <A point on the normal to the middle surface will have a 
constant distance Z from this surface, both before and 
after deformation. 

Accordingly, 


Ll wen. Z ow. WwW. 
Geni ON ee. 
e A d¢ ~ MATZ) Nae AtZ 4 


Substituting Equation (3) into the definitions for N and 


M, Equations (2a) and (2b), leads to 





a etm % 


Ws w| (5) 


K 


me 


Qy 





M = ye 
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ei 
where D=Eh eon SL 


The linear equilibrium equations for the element follow directly 


from an examination of Figure 3. 











Theseus 

Se — 2 arene Pah . (6a) 
So atl” oT a! 
Q _< q (6c) 


where rotary inertia is neglected. 


Eliminating the radial shear Q from these equations leads to 











ON ~ as OM. = eAh ay 
ys A 9S¢ (7a) 
2 
| oM 2 
2 zr = PAh 2 (7b) 
Set N= Cdn eM 
Substituting Equations (4) and (5), relating the stressresultants 
to the displacements W and V, into Equations (7a) and (7b), gives 
D/XV _ aw\. o°V 
ASH 39) 7 Ah Sy (8a) 





(SK +23 34 +w)- W367 PPP Ab ae = oa 


These are the differential equations that describe the free 


vibrations of the thin, two-dimensional, elastic ring. 
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CHAPTER IIT 


SOLUTION FOR THE IMPULSIVELY LOADED RING 


In order to adequately test the proposed method for increasing 
the critical time step, an impulsive load that causes both large 
bending and membrane stresses should be employed. Consequently, 
the impulsive load illustrated in Figure 4 is taken as the loading 
on the ring at time © = O.m@@The impulse is, applied to the top half 
of the ring and is characterized by a cosine distribution directed 
radially inward. 


The initial conditions for the radial cosine impulse are given 


by 
W(¢,0)=0 V(b,0)=0 
ot, = 0 
SW (¢,0) = -v, Cos > EeosZ 
= 0 -ELeesF (9) 
where ve tS 


eh 


Since the ring is unrestrained, a net momentum will be imparted 
by the radial velocity acting on the upper half of the ring. This 


will cause a rigid body translation of the ring, and hence, unbounded 
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FIG. 4 
IMPULSE LOADING ON THE RING 
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values for the displacements V and W. To avoid this situation, the 
initial conditions must be modified in such a way that the rigid 
body motion is eliminated, while preserving the internal effect of 
the impulse. Reference 2 defines the following modified set of 


initial conditions for this loading 


W(d, 0)=0 Vio, 0)=0 





ove _ Ve 
7 ($,0) 4 SIN. ® 


ow ie 3m 
(>) 0)= 3h cos > HET cokes 3: 





# COS -Bepe F (10) 


With this set of initial conditions, the ring will remain fixed in 


space; and hence, the displacements V and W will measure actual 


deformations of the ring. 


The following solutions for W and V are assumed 


Wl, T ) = > Ay (T) Cos np (lla) 


V(>,T) =) Br (T) SIN nd (11b) 
n=! 


. * 
(WyT 


LG), T 
B,(T)= b,€ (12) 


Ay(T)= &n e 


*a and b= are arbitrary constants. 
n n 
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Based upon these solutions and the given initial conditions the 


expressions for the bending moment and normal force are* 





2 = Wri [2 2 
M, = -mM Awe = SIN®T + Bil) m=n SIN £,W.T 
N=24,... 
2 .¢ 
‘ol SIN gnW.T| COS Nq (13a) 


No= NA@T s_s~at + —% SINVZ WT COS & 





DY, 22 

= MW + | 
2 (a1 ee 
ee a | ae 
(n8-) (48a 98) £, SIN f, 

224... 

+ 9n SIN 9n &7| Cos nd - aM, (13b) 
where ot = b MK 8 d), _ | D 


A Veh 


M, and Ny are non-dimensional representations of the bending moment 


© 


and the normal force, both of which may be solved without specifying 


the impulse magnitude, I The other terms will be described later. 


= 
Equations (13a) and(13b) characterize the dynamic response of the 


elastic ring due to the specified impulse. Before proceeding further, 


however, a physical interpretation will be giveto these expressions. 


*These expressions are derived in Appendix A. Since this is 
the original ring, § = 1. 
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Consider, first, the various mode shapes excited by the impulsive 
loading. Figure 5 illustrates the lower modes of vibration for the 
ring. The lowest mode, n * 0, is characterized by a radial displace- 
ment only, where Wis the same for all points. A small change tn 
curvature in addition to the middle surface straining will result; 
and hence, N, and M, will both have contributions from this mode as 
can be seen from Equations (4) and (5). The next mode, n = lL, 
involves middle surface straining only. The shape of the ring 
remains the same, with no change in the effective curvature of the 
ring. Hence, this mode will contribute only to the normal force, 

No. For n> 1, only the even mode shapes appear in the expressions 
for the stress-resultants. This situation arises from the modified 
initial conditions imposed on the ring. At time t = 0, positive 
radial velocities are required at the top and bottom of the ring, 

@ = Tr and 0, while zero radial velocities are prescribed at the 
sides, @= 22, 


the even modes can account for these velocities and hence, are 


Referring to Figure 5, it can be seen that only 


the only ones excited by the impulse. 

Finally, examination of Equations (13a) and (13b) reveals that 
two frequencies, f We and g, We» are associated with each mode 
of vibration of the ring. These frequencies arise from the fact that 
the ring element, dad , may displace independently in two separate 
directions. From Appendix A, the frequency coeffictents, ee and 


B? are approximately given by 
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LOWER MODE SHAPES 





f, = n(n?-1) Er aag 
Gn = J n?t] 


Since a thin ring is assumed, “is very small, and ome go The 


cle 


frequency coefficients, Sa, and g > characterize, respectively, the 
radial and tangential vibrations of the ring. Due to the relative 
magnitudes of fT and 8 the expressions for M, and N, will each be 
Gom@nated by one Of sGhe two sine terms %of the series. @Hence, M. 
will vary at a frequency given approximately by £ We ; while N, 
will be characterized by a much larger frequency 8 Wo - 

To obtain numerical data from Equations (13a) and (13b), it is 
necessary to select physical data describing the ring. Figure 6 
outlines the physical characteristics for a thin ring composed of 
aluminum. The frequency coefficients, rT and go have been calcu- 
lated on a digital computer at 20 micro-second intervals over a 
periodmot time weorresponding to one complete cycle for M, . Streseae 
on the inner and outer edges of the ring were also calculated from 


the following non-dimensional relationships 


(outér edge) S. = —- Oo, = N, sm > M, (15a) 
(inner edge) G) = A Woe TT sr =N - ate M, (15b) 


E V, 
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16980 rod./sec. K = 1.43 x10° Ibf.-in. 


fn = 
D = 2.74 x 10% Ibf./in. p= 3.14 x10" eae 
in. 
E = 10.9 x 10° Ibf /in? Q:= 354x10° 
FIG. 6 


PHYSICAL RING DATA 


aah 


‘ ' : 
The maximum values for Me 2. News and OG, found over the time sampled 


have been listed in Table I.* 


Figures 
ring through 
Table I that 
and has been 
is identical 


The symmetry 


7 through 10 illustrate the dynamic response of the 

Eremtirst half cycle of M, .** It,Should be noted from 

the response at <>) = tL 2 is of secondary importance 
4. “4 

omitted from the curves. Also, the response at @= 77 


in character to that at @= QO and has been omitted. 


of the even modes make it unnecessary to determine 


solutions for ring locations in the range 77 ¢€ @ €< 277. 


ae ( ane 
*The maximum @, was observed to be within a few per cent of 
the maximum on found. 


*kA positive M, indicates a negative bending moment. 


28 


TABLE I 


MAXIMUM VALUES FOR THE STRESS-RESULTANTS 





] 

P No nn Mr oe 
0 202 “ome BA 
us 1.80 26.5 1.99 
A 

us eet. 125.8 2.57 
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oT ie 26.5 1.80 
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CHAPTER IV 
SOLUTIONS FOR THE IMPULSIVELY LOADED 


MODIFIED RING 


The assumption that the ring's inertia characteristics are 
variable with regard to tangential vibrations will now be considered. 
The ring equations must be modified to account for this assumption, 
and modified solutions analagous to the original ring are obtained 
in Appendix A for the impulse loading described in Chapter Ii. 

Figure 11 illustrates the assumed inertia forces associated 
with the ring element, dd - A Simple adjustment of the ring 
equations will account for this situation. The following expres- 
sions for the stress-resultants, derived in Appendix A, are obtained 


for the modified ring 


— B+! 
7 - eur 
Ms, = SIN GT + > 2a |e SIN 4, 03,7 
8(F?-34) - 97 ) £, ’ 


N= 2,4,--- 


2 —Z 
vg eel S] 


a N Qn Wer | COS NP vite 
n 


MN TT : 
N= -SIN@T + wireval SIN (i+ % 7 COS $ 


[iG 


= Yt 


ya Cs)) | 
1 OIE Oy) | 





—f Sif, WT 


= 
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34 
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FIG. Il 
ASSUMED INERTIA FORCES 


Se 


wa 


To determine the effect of the factor 8 on M, and N, , consider 
first the change in the frequency coefficients, i and go: As deter- 


mined in Appendix A, these coefficients are 


— 2 (am 4 
= a) oe (18) 
fi, (n?-1) iPep : 
= iit n*+yx 
= \ (19) 


For this analysis, it will be assumed that ¥ takes on values between 





l and 4. Therefore, i will not be significantly affected by 6 , 
especially for the higher mode shapes. However, em will be very 
sensitive to a change in 6 , as can be seen in Equation (19). Figures 
12 and 13 graphically illustrate the error in the modified frequency 
coefficients as a function of 8. A physical interpretation of the 
factor § will verify the above results. As 8 increases, a progres- 
sively heavier system is seen in the tangential direction, giving 
rise to a response of lower frequency, 8 GW, » The radial vibrations, 
characterized by ie GD,» should not be significantly altered by ¥ 
since the density in this direction is always given by the original 
ring density, ‘2 ‘ 

Values for the modified frequency coefficients have been calcu- 
lated and listedsin Appendix B for the ring defiméd in Figure 6. 
Again, it can be seen that eet Lil a M. and N, , then, will be 
approximately described by sinusoidals at frequencies E We and B Wo » 


respectively. 
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FIG. 12 
CHANGE IN f, WITH ¥ 
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FIG. 13 
CHANGE IN §, WITH Y 
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Examination of the expressions for My and N, reveals the 
dominant terms. For the bending moment, Equation (16), the lower 
frequency sine term will dominate. The magnitude of the coeffi- 
cient for this term will decrease as the product, 9& a increases, 
corresponding to increasing values for ¥ . Thus, the time 
response for the bending moment will decrease in magnitude as YX 
increases. The wave form for the response, however, will not change 
significantly. Examining Equation (17) in a similar manner, it can 
be seen that both the wave form and the magnitude of the time response 
for No will change as ¥ increases. 

Figures 14 through 19 illustrate the curves for M, , No, and 
ore for three values of § at p *Q . Figures 20 through 22 show 
M, for @ = ue .* The values obtained for M, and N, at @ = - and 
aa) = FP were found again to be small compared with other Locations 
on the ring and have been omitted. Also, the curves for > = TT 
have not been included because of the similarity with the response 
ror d = OF 

Tables Ii through LIV list the maximum values observed for the 


Z ee | 
solutions of M, , Nj and O, . The corresponding values for the 


unmodified ring have also been included for comparison. 


*N, and G,/ at p* 4a are not included because of their small 
magnitudes. 
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TABLE [II 


MAXIMUM BENDING MOMENTS 


=| S=2 y=3 Oz=4 

-129.6 -129.5 -124.9 -117.7 

1256 Le..3 L1Omes TOS a2 

-129.6 -129.5 -124.9 -j1i7./ 
TABLE ill 


MAXIMUM NORMAL FORCE (MAGNITUDES) 


Se Y=2 §e3 §=4 

ape AE 230 Ze a0 2.48 

1 1.48 1.14 Laie 

es. Oe 2.43 2.46 Z.48 
TABLE IV 


MAXIMUM OUTER EDGE STRESS 


X=] = 2 x3 Sud 
-3.45 -3.49 -3.28 =32.5/ 
Pb OM PARES \ 2 a0 Zale 
ao. 2 / -3.09 -3-31 -2.98 
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CHAPTER V 
CONCLUSIONS 

A comparison of the response of the original ring with the 
results for the modified ring with XY = 2, 3, and 4 can be made using 
the data presented in Tables I through IV and the response curves 
shown in Figures 14 through 22. A study of the time response for 
the bending moment, M, and My ; reveals that the solutions to the 
modified ring, in general, range from fair ( § = 2) to poor ( gf = 4) 
approximations to the solution for the original ring. Average errors 
range from 10 per cent for = 2, to about 40 per cent, for § = 4. 
However, the percentage error associated with each modified curve 
is nearly constant, and the wave form is identical to the unmodified 
curve. Jt would appear, then, that the data from a modified solu- 
tion could be increased by some constant factor to yield a very 
good approximation over the first quarter cycle. 

Based upon the above comparison, the proposed method appears to 
introduce unacceptable errors into the solution for ¥ > 2. However, 
the actual time response for the bending moment is of little signifi- 
cance in a dynamic analysis, since the maximum value of bending 
moment is usually the desired result. The maximum bending moments 
observed for all solutions have been listed in Table II. Based on 
these values the solutions for the modified ring are very good 
approximations to the results for the unmodified ring. Errors in the 
maximum bending moment of about 4 per cent and 8 per cent can be 


seen for § values of 3 and 4, with essentially no error observed 


at 1G eae 
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The effect of a heavier system in the tangential direction 
is readily observed in an examination of the various normal force 
response curves. Both the wave forms and the maximum values are 
affected by a change in the factor 6. As anticipated, the tangen- 
tial vibrations decrease in frequency as 8 increases. This can be 
seen from the spacing of the beats. Despite the obvious changes 
in the normal force curves for the modified ring, the maximum peak 
values, listed in Table III, are not significantly altered as @ 
changes. 

The primary concern of any dynamic analysis is usually the 
maximum stress that occurs in the ring. Since both the bending 
moment and the normal force contribute to the stresses, the 
maximum value will be found on one edge of the ring. The relative 
importance Gf ~eaen coOntribwtor, M, and N, » can readily be s@en 
from the stress curves found in Figures 14 through 19. Table IV 
lists the maximum non-dimensional values found from these curves. 
The bottom of the ring, d =QO , is seen to be the critical stress 
location. A non-dimensional value of -3.45 is the maximum stress 
calculated from the solution for the original ring. In comparison, 
the modified ring solutions have values of -3.49, -3.28, and -3.57, 
all of which are within 5 per cent of the original ring result. 
Consequently, the assumption of an artificial tangential density 
has a negligible effect on the ring's maximum dynamic response for 


the impulsive loading condition considered. 


B18 


In conclusion, the proposed method for reducing the stability 
restrictions presented in Chapter [ appears to be valid. The same 
method when applied to the finite difference approach should produce 
results very similar to those presented here. Thus, because of the 
assumed artificial tangential inertia, an accurate finite difference 
solution should be obtained with the computation time reduced by 
about 30 per cent for § = 2, 40 per cent for v= 3, and 50 per cent 


for §= 4 from that of the unmodified ring analysis. 


ay 
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APPENDIX A 
SOLUTION FOR THE SPECIFIC IMPULSE LOADING 

The derivation for the solution to the radial cosine impulse 
is outlined in this appendix. It should be noted that the develop- 
ment is identical to that of the exact solution when the factor 
eal. 

Assuming that the density in the tangential direction is larger 
by a factor § than that associated with the radial direction, the 


following modification in the ring equations will result 


ae - 4). se An SY, (A-1) 
V 


ce a ee 


Solutions for W and V were assumed earlier in Equations (lla) 





and (11lb) and are repeated below. 


w(.T) = ye A, (T) Cos n@ 


(a3) 
v( ,T)= 3 B,(T) SIN NP 


The symmetric and yee forms for W and V arise from the 
symmetry of the impulse and the definition of the V, W coordinate 


system. 


oe) 


Substituting the assumed solutions into the modified ring 
equations, Equations (A-1) and (A-2), and setting the two algebraic 


expressions equal to zero for each value of n lead to 





Dn( Ay - By) A? 2 8n = (A-4) 
a n mes C oT? =O 

2 2 2 

D( An - By) + An ae (N*~1) +heh or =O (A-5) 

ror no) | ae , oer OS 
Harmonic forms are assumed for An and BL 
Ay (T) = on ial 
(A-6) 


‘y (OnT 
By 6p = by € 
where a and be are arbitrary constants. Substitution of these 


assumed solutions for the time functions into Equations (A-4) and 


(A-5) yields the following two expressions from which the common 


time term, of Ont , has been eliminated 
Dn (a, -nb,) + sMeh wy by = 0 (A-7) 
2 
D(a, -Nby ) ie Ay + (n?-1) -Neh Wy Oy, =O (A-8) 


The coefficients a, and b cannot be determined from the above 
equations since a third unknown, G)y, is involved, However, they 


may be expressed {n terms of one another and eliminated from hoth 


expressions. 


From Equation (A-7) 


_ heh wn 














Qn = b N (A-9) 
n 
Dn 

Replacing a in Equation (A-8) with the expression above, the 

n 
constant ba is common to all terms and may be factored out. A 
single equation in the unknown wo, results 

4 2 2 | 2. 
8 | Dn x xa (n°-1) @,, 2 ,\4 
— | — |} — | Ft NI | vt —NAN -/ = -10 
NY \W, n n Wo ( ) 0 Ca 
where 





: ae a 
OG), rN Ph of DAZ 


The modal frequency wo for the time functions AL and BT may 
now be determined. Although four values for GQ) will satisfy 
Equation (A-10), only two, the positive roots, have physical signifi- 
cance. These roots indicate that each mode of vibration is character- 


ized by two distinct frequencies 


(ASI) 


* The bars indicate that the frequency coefficients are 
functions of @. 
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where 





Z ; [(n%+5) +ocy (n?-1) -\/(ntes) 20 (n2-1) (sn@-52)+ v7? (ne) i 
n= RA ANTES 
20 





(A-12) 
ge lentvs) +c (eft) + Vinee -2e (PF (un= 8) Foe (I | 
n ie 


The frequency coefficients A and 8, may be simplified in Light of 
the thin ring assumption. For small values of &*,~ these 


coefficients are approximated by 


Pl as 2 ot 
fy = N(N=1)\s 





(A-=13) 


nee 


ae 
The initial values for the time functions must now be determined 
in order to solve for the constants a, and bi At time T # Q, the 
expressions for the displacements and velocities are siinple Fourler 
series expansions in the single varlable d. The Fourler coeffi- 
clents for these series may be evaluated to yield the following 
initial conditions for ae and BCT) 


A (0) 2 () id 30min 2 pees 


“For the problem considered, K* 3.54 x 10’. 





ade 

abo) = -$ 
0 

By (0) =0 

SS (0)= = 

$80 (0) = O 


a 
I) Nigel Ce 
one 
na 3, Ss 
ie ol ee 
mee, 2, 3, 4, 


(A-14) 


a 


The time functions may now be determined from the above initial 
values. For the cases n = O and n = 1, a simple substitution into 


the modified ring equations yields 


ma. SIN WeT (A-15) 
77 W 


-) 


A, (T) 


and 


-~V, 
4 [+k We 


A,(T) SIN I+ 4 WeT 


(A- 16} 


Vo 
B,(T) a 4)i+h we SIN I+ G),T 


The solution for the higher mode shapes, n > 1, requires further 
explanation. Recalling that the modal frequency Wy was character- 
ized by four distinct roots, there are four separate solutions for 
AL and BY which satisfy the ring equations. The sums of these 


solutions are also solutions. Therefore, the general time functions 


for n > l are 


Fat — -FrwT — jGnw,7T — -iGnwr 
tts ote oe 


Ay (T)= On nz © + ay3z € nd. 


(A-17) 
Ts en Gn WoT =? Bn WT 
By (T)= by, € ° + bye € + by, e 8 "+ Png © . 
Although eight constants appear in the above equations, only four 
are independent. As determined earlier in Equation (A-9), the b 


constants may be replaced with expressions involving the a constants. 


Applying Equations (A-13), these relationships are 
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Y) — 


aoa oer 
ni On n7- yf? ) by2 = a yy2_ yf? 


o | 
i 
| 


(A-18) 


o | 
= 


ne * Ong oer ae] i Bnet Snel ees 
Substituting Equation (A-18) into Equation (A-17), the expressions 
for A_(T) and BAT) contain four unknown constants to which four 
initial conditions may be imposed. 

To apply the initial conditions, Equations (A-14), the time 
functions must be expressed in terms of sines and cosines. Making 
the appropriate expansions for the exponentials and regrouping the 


@ constants into new constants, 4a, Equations (A-17) become 
An (T) = Gy, cos fy) WoT + Gy2 SIN fy, WT 


B,(T) = an COS FOLT + Gwe SIN ty 07 





ste i G.. COs Qn&.T + Qi, SIN Gn OT 
Vy- 2 n3 Ni4- 
-¥9n 

A solution may now be determined for the 4 constants. Applying the 
initial conditions, the cosine terms must vanish. Thus, 

== = = Bey 

Qn, i. Qn 7 bn, i. by3 = & (A-20) 
The coefficients for the sine terms are given by the following 


expressions 


a2 


ov, 2" (8h? - ni?) 


rol 
ie 8 W,17 (N?-1) (F,7 - an) F, 
1 = eee. 6: 
"n+ = 
= 2a (i)! (rn? we) aa 
Mee ¥ Wo, TT (1n*-1)(F2- a2) a, 
aie = (ely) = O a= 3, le 
and hence n 
= ay, one = 
re YO. = NER 38 ) Fy die. 
(A=2) 
= ie 2 Ve (- yet! yn 
Pg 6 Gd. 1 (ni? -1)(F,2- FF) Fn 
o> Wem = O n> 3, 2, 
(A-24) 


The time functions may now be written directly as 


+ | 
7) aigGne a ayer -n* n?-v5n 
A,,(T) ¥oa.m (n@1)(GE2 38) a Sut 2 coer @ Pe = SIN Gn oT 


i. =ecy F 
n(T)= ya, 71(n?-I) (2-3 2) | p SIN pigoel <t aet a, SIN Gn We 7 
even modes > 2 (A-25) 
Ay (T) = By (T) = 0 odd modes 33 (A-26) 
These expressions may be substituted into Equation (A-3) to obtain 
the middle surface displacements, W and V, as functions of time and 


ring position. As determined earlier in Chapter II, Equations (4) 


and (5), the stress-resultants M and N are given in terms of the 


displacements W and V by 


N = 2 (SY —w) - eee + w | 





ays Me \ Spe 
(A-27) 
_ Ky] ew 
M mabe + Ww] 


Thus, the solutions for W, V, M, and N due to the radial cosine 


impulse become 





W he =o NaGae — 7h Feral INV I+ wr COS p 


O© 


Worl r= 
+) oo a 7 i Saf, Qe 
n 





¥ (y2-] 2 
= 24... n Ny - 3 tr 
7 | 
ues: J SIN ce nod (A-28) 
Y) 
Oo = —_i__aeGiiyi ieee ana 
V V, 4 1+! ft Pp 
eS +1 
2 (-1)” 
+ eo SIN £, WeT 
) ¥(n?-1)(£7-3 a|e 
N=24,--- 
es 
ee SIN Sn a SIN Np an 
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rv, A’ ar w) 
NAL a | KV, 
ee pe ly oe a 
eS Co.7 +) oe Maat SIN 4, WoT 
geeaa. Se - 5) n 
» W= 2 ow swt (A-30 
%, Sn WT} COS Np ) 
N, = NAMw& . —SiIN @iy + get SIN + CRT COS 
DV, ee 7" I+% f 
oO + | 
\ al S46) 5 i) a ae i 
2D We (n?-1(F2 a8) -£, SIN ty WoT 
+ an SIN On sar |cos Nb _« M, (A-31) 
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APPENDIX B 


CALCULATED FREQUENCY COEFFICIENTS 


TABLE V 


FREQUENCY COEFFICIENTS ( & = 1) 


0.016 
0.088 
0.208 
UIA a7 ha 
Umea 
Oe S57 
1.169 
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TABLE VI 


FREQUENCY COEFFICIENTS 


i 
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FREQUENCY COEFFICIENTS ( 6 = 3) 


TABLE VII 
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370 
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FREQUENCY COEFFICIENTS ( § = 4) 
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O13 
081 
moo 
PoOy 
200 
848 


161 


By 

414 
2:36 
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123 
099 
082 
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